Critical quench dynamics in confined systems 
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We analyze the coherent quantum evolution of a many-particle system after slowly sweeping a 
power-law confining potential. The amplitude of the confining potential is varied in time along a 
power-law ramp such that the many-particle system finally reaches or crosses a critical point. Under 
this protocol we derive general scaling laws for the density of excitations created during the non- 
adiabatic sweep of the confining potential. It is found that the mean excitation density follows an 
algebraic law as a function of the sweeping rate with an exponent that depends on the space-time 
properties of the potential. We confirm our scaling laws by first order adiabatic calculation and 
exact results on the Ising quantum chain with a varying transverse field. 



Recent beautiful experiments with ultracold atomic 
gases have revitalized the studies of nonequilibrium as- 
pects in strongly correlated quantum systems pQ. The 
main peculiarities of the dynamics of ultracold atoms are 
low dissipation rate and phase coherence over very long 
times [TJ such that the dynamics is very well described by 
the usual quantum unitary evolution of closed systems. 
An important problem that has received much attention 
is the situation where the parameters of the quantum 
many-body system are varied in time such that it reaches 
or crosses a quantum critical point 2J . The divergence of 
the intrinsic relaxation time close to the critical point in- 
duces non-adiabatic evolution no matter how slow is the 
rate of change of the Hamiltonian. If the system is ini- 
tially in its ground state non-adiabatic transitions toward 
excited states lead to the presence of topological defects 
in the final state [3H7]. For example, driving a quantum 
system from a paramagnetic to a ferromagnetic phase 
through a critical point generates a final state given by 
a superposition over excited states carrying finite ferro- 
magnetic domains separated by kinks or domain walls. 
For a slow driving rate the density of defects is a univer- 
sal scaling function of the driving rate as in the classical 
Kibble-Zurek (KZ) mechanism [5]. This may be of impor- 
tance in the context of adiabatic quantum computation 
[§] where adiabatic evolution is proposed to transfer the 
system from an initial state to a computational nontriv- 
ial state. If one is forced to cross a critical point in order 
to generate the nontrivial state, inevitably the crossing 
will result in the generation of excitations (defects). The 
optimal time ramp needed to drive the system through 
the critical point has to balance the unavoidable genera- 
tion of defects and the time needed to cross the critical 
point |7]. 

In many real physical situations the critical systems 
under study are immersed into an inhomogeneous field in 
one or several directions (for example the gravity field) 
whose main effect is to smooth out the critical singular- 
ities [T0HI2] . This is particularly relevant in the context 
of ultra-cold atoms where parabolic trapping potentials 
are used to confine the atomic cloud. Removing or load- 
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FIG. 1: Sketch of the slow sweep of a power law confin- 
ing potential. The state of the system (red) deviates from 
the adiabatic ground state (blue) as we approach the critical 
point associated to a vanishing potential. The small stripes 
represent the spectrum at different times. 



ing smoothly in time such a power-law trap will lead to 
a final state carrying a nontrivial density of defects, de- 
pending on the shape of the trap, as soon as we are get- 
ting closer to a critical point. In this letter we analyze 
the coherent generation of defects in many-body quan- 
tum systems after slowly sweeping a power-law confining 
potential close to a critical point, see figure[Tj We develop 
a general scaling scenario from which we derive the scal- 
ing behavior of local (such as the local energy density 
or order parameter) and global quantities (such as the 
density of defects). We derive in particular the optimal 
non-linear time ramp minimizing the generation of de- 
fects for a given total sweeping time. The scaling theory 
is tested on the Ising quantum chain with an inhomo- 
geneous transverse field playing the role of the confining 
potential. 

Scaling theory Consider a system which has at 
zero temperature a quantum critical point, controlled by 
a scalar field h, separating a symmetric phase from a 
broken-symmetric one. Assume that close enough to the 
critical point the control parameter h deviates in one di- 
rection from its critical value h c with a power law 



S(x,t) = h(x,t) - h c ~ g{t)\x\" 



(1) 



The amplitude g(t) of the deviation is driven externally 
and varies in time as g(t) — v |£| Q sgn(£) where, without 



loss of generality, the rate amplitude v is positive. The 
space and time exponents ui and a are both positive. 
Notice that within this parametrization the quench dy- 
namics connects the two phases (d > and 5 < 0) by 
crossing the homogeneous critical point at t = 0. 

The inhomogeneous perturbation (fTl) introduces a 
crossover region in space-time {x, t) around the critical 
locus (0, 0) with characteristic length-scale £ and time- 
scale r. Indeed, starting in the far past in the ground 
state |0(g(— oo))) the system evolves adiabatically toward 
the instantaneous ground state \0(g(t))) as far as it is 
protected by a large energy gap A(t). In the instanta- 
neous state \0(g(t))), the spacial power-law variation of 
the control parameter S introduces a typical length scale 
£{t) around the spatial critical locus (here at X = 0) [TT] , 
The typical length can be obtained self-consistcntly from 
the correlation length relation £(t) ~ 5(£, t)~ u . One finds 
£(t) ~ \g(t)\~ v / ( 1+vu > diverging for a vanishing devia- 
tion amplitude g. Here g(t) plays the role of an effec- 
tive deviation to the critical point with effective expo- 
nent v g = l/i/g = ^/(l + vlo). As time runs toward zero, 
the gap closes and the relaxation time ex 1/A grows up 
to the point where the adiabatic approximation breaks 
down completely. Getting closer and closer to the criti- 
cal point the dynamics switches to a sudden regime. For 
larger times, after the critical point has been crossed, 
one recovers again the nearly adiabatic regime. This 
is the Kibble-Zurek (KZ) scenario. The typical time- 
scale r around the critical locus where this break-down 
occurs can be deduced from the self-consistent relation 
t ~ £{t) z , with z the dynamical exponent. It leads to 
T r^ £ z r^ v - z /yv wn ere y v = y g + za = (1 + v(uj + za))/v 
is the RG dimension of the perturbation field, such that 
under rescaling by a factor b the amplitude v transforms 
as v' = b Vv v. Since to and a are positive, the perturbation 
is always relevant (y v > 0). 

Local and bulk quantities The profile of a local quan- 
tity tp(x,t, v), like the local order parameter or the en- 
ergy density, with scaling dimension x v transforms under 
rescaling as f{x,t, v) = b~ x ^(p(xb~ 1 ,tb~ z ,vb Vv ). Taking 



l/Vv (x £ (x r 1 / 2 , one obtains 

<p(x,t,v) = v x ^ /v ^(xv 1/y ",tv z/Vv ) 



(2) 



The prefactor exhibits the trap-size scaling ip ~ l~ x f 
associated to a finite size system with £ ~ v~ 1 ^ Vv [TT]. 
The unknown scaling function $ encodes the propagation 
in space of the critical phase (close to x = 0) into the 
initial one. See [TJ5] for detailed studies on that issue. A 
bulk quantity or a spatial averaged profile Tp(t, v) satisfies 



Xp(t,v) = v x * /y "$ (tv z/y A 



x <p/ z (f) 



(3) 



As an example, the spatial averaged energy density, with 
scaling dimension x e — d + z, should behave after the 



Density of defects An estimate of the density of defects 
n generated by crossing the critical point with the non- 
linear ramp g(t) is obtained from the length scale identi- 
fied above through the relation n ~ £~ d ~ v d ' Vv . To be 
more precise, equating the relaxation time To/A(£) with 
the typical time-scale A(i)/|A(i)| at which the Hamil- 
tonian is varied and assuming that the gap scales as 
A(i) ~ fto|# /y9 one finds for the typical (Kibble-Zurek) 

/ \Vg/Vv _ . 

time-scale tkz ~ o — v z ' Vv . From the relation 

n ~ [A(T KZ )] d ^ z we get the behaviour 



)da/y v 



4/vv _ 



(zj6) 






(4) 



with 6 = ^w 1/q ~ 1/T (T defining the temporal win- 
dow of the quench protocol) and 7 = av g = a/y g . As 
expected the density of defects n is smaller for larger val- 
ues of the protocol time window T necessary to reach a 
final value gf from the initial go one. If one wants to 
minimize the generation of defects, the switching of the 
trap should be as slow as possible. However this can lead 
to extremely long protocol times T and be counterpro- 
ductive, as for example in quantum computational issues 
where one looks for a compromise between the produc- 
tion of excited states and short computational times. To 
achieve this compromise one may look for the optimal 
power-law time ramp protocol that minimizes the defect 
density n at a fixed duration T. Extremizing J4J with 
respect to 7 = a/y g for a given S ~ 1/T one finds 



lo P t = -W ( - x 
z \ eo 



(5) 



quench to the critical point as e 



v 



(d+z)/y„ 



where W(x) is the Lambert W function defined through 
x = f{W) = We w . For a given trap shape (space expo- 
nent uj fixed) the optimum time exponent takes the value 
a op t = lopt/v g = 7opt(l + vu)/v. The result of [7 is re- 
covered with u) = by expanding the Lambert function. 
We see that loading a trap potential changes significantly 
the homogeneous optimal protocol by increasing the ex- 
ponent a by a factor (1 + i/cj) (close to the critical point 
one has to be slower than in the homogeneous case). 

Dynamics in the Ising chain The Ising quantum 
chain is a standard theoretical laboratory for issues re- 
lated to quantum phase transitions [2] and it has been 
already studied in various sudden quench dynamics con- 
texts [14]. In this model, an inhomogeneity analogous to 
the one arising from a trapping potential in particle sys- 
tems can be achieved by considering a time-dependent 
inhomogeneous transverse field. The Hamiltonian is 

L-l L 

n(g(t)) = -2 E CT >-+i " 2 E &»GK*)K> ( 6 ) 

n— 1 n— 1 

where h n (g(t)) = 1 + g{t)n u , with g(t) = v\t\ a sgn(t). 
Notice that the spatial critical locus has been set at the 



left boundary of the chain. For the critical unperturbed 
system, g = 0, the dynamical exponent z — \ and the 
correlation length exponent v = 1. 

If the rate of change of the Hamiltonian is slow enough 
one can use a nearly adiabatic approximation in order 
to describe the actual state |*(i)) generated from the 
initial ground state |0(io)) by the unitary time evolu- 
tion. Introducing the instantaneous eigenbasis {\k(t)}}, 
such that H(g(t))\k(t)) = E k (t)\k(t)) , one obtains from 
standard perturbation theory the "one-jump" expansion 
|*(i)) ~ e- id °(«°-'>(|0(i)) +E k ^o a k(t Q ,t)\k(t))), where 

$k{to,t) = L ds E}~ is the dynamical phase. The transi- 
tion amplitudes ak(to,t) are given by 



a,k(t ,t) 



9(t) 



9(*o) 



8u>ko(g) 



(7) 



with SuJko(t) — Ek(t) — Eo(t) the Bohr frequencies and 
$ko = flk — fio a dynamical phase factor. Under a Jordan- 
Wigner fermionization followed by a canonical Bogoli- 
ubov mapping the Ising Hamiltonian (16J) is diagonalized 
in terms of the Dirac fermionic algebra {gt(g) , Vq(g)} = 

S pq and reads H(g(t)) = J2 P e p(ff) [^(flO^pCff) ~ V 2 ] 
where e p {g) are the instantaneous one-particle excitation 
energies. The perturbation d g T-L{g) in Q takes the form 

dMg) = lY. X p^M(9)+Vp(9)M(g)-v q (g)l (8) 

p,q 

with X% q {g) = Y^n^pi^d^ipqi^g) where the (j> and 
i\) are the Bogoliubov coefficients. This perturbation 
induces transitions from the ground state to the two- 
particles states \pq(g)) = »?J (5)^(5) 1 0(5)) only These 
fermionic particles are the basic topological defects in 
the Ising chain. At time t, the density of defects in 
the mode q is given by the instantaneous population 

n q (t) = (^(t)\vl(g)v q (g)\^(t)) ^ ^E P K q (t ,t)\ 2 . 

Linear spatial perturbation For a linear spatial modula- 
tion, that is at u — 1, we have an explicit solution in the 
thermodynamical limit L — > 00: the Bogoliubov coeffi- 
cients are given by the wave-functions \p ( U P to normal- 
ization) of a simple harmonic oscillator |12j . The matrix 
elements (pq(g)\d g 'H(g)\0(g)) in (IT]) are then proportional 
to the position matrix elements of the harmonic oscilla- 
tor {xpi uxq), such that the only non vanishing transition 
amplitudes a pq (to, t) are those with p = q ± 1. 

Plugging the exact solution into ^ one obtains a 
closed expression for the amplitudes a pq . However, con- 
trary to the spatial homogeneous case (to — 0) where the 
integral ((7| converges at the critical value g = 0, here the 
linear spatial inhomogeneity modifies the dependence on 
g of the integrant to a g^ 1 behavior leading to a logarith- 
mic divergence at g = 0. This divergence is caused by 
the square root dependence on \g\ of the excitation spec- 
trum e p = \g\ 1/2 J\p + 1 + sign(g) with p = 0,1,2,... 
[12) . Consequently, the first order adiabatic expansion 



breaks down at the critical point 5 = (at t = 0). Nev- 
ertheless, for quenches that do not cross the critical point 
(the starting and the ending point of g are on the same 
side of the critical locus) one can still use |7| giving 



lo-pqih,^ 



Ui 



>4fl 



A Pp (\go\, \g(t)\) 



(9) 



/«/ 



with G pq = J2{p + q) + l + sgn(g)[<5 pg _ 1 - 6 pq+1 ] 



ttpq = \9\ 1/2 ^pqfii.g) and Ppq = -20; 



The function 



PI a+2 



sgn(ff). 



LCi r / 2 + a \ / 2 + a \ ' 

A p (x,y) = j^ [E a (ipx^j - E 1 (ipy^)\ (10) 

is expressed in terms of the exponential integral E^z) = 
J^°dtt- 1 e- t for |Arg(z)| < tt. 

In order to compare the analytical result (p)]) with 
our scaling predictions consider first the case where the 
quench starts far away from the critical locus, \to\ 3> 1. 

In that case Ei [ip\go\^^ ) — and the function A(x,y) 

is dominated by Ei(ip\g(t)\^^) leading to the expected 
scaling behavior n q (t,v) = f(v\t\ 2+a ) with y v = 2 + a, 
see (pi). On the contrary if the initial trap amplitude 
go is small enough one observes the non-homogeneous 
behavior n q (g ;t,v) ~ fo(go,v) + fi{\t\v 1/v ")- The ex- 
pected scaling behavior pi is broken by the presence of 
the boundary term /o which accounts for the high correla- 
tions in the initial ground state \0(go)} (since for \g \ <C 1 
the system is nearly critical). 

For a quench at or crossing the critical point the sit- 
uation is more complicated since, as stated before, the 
(un-normalised) perturbation formula nty leads to a di- 
vergence at t = 0. However, for a finite size chain the 
energy gap Suika stays finite at the critical point which 
wash out the critical divergences and one can perform a 
finite size scaling study. 

Finite size scaling analysis for general to We start far 
away in the disordered phase {g — 1) and drive the sys- 
tem to the critical point (g = 0). The density of defects 
is obtained numerically by exact diagonalisation of finite 
chains with up to 256 spins. The finite size data are then 
extrapolated to the infinite size limit. The results are 
reported on figure [2] for different spatial and temporal 
exponents w and a. First of all we observe that the finite 
size values are always smaller than the extrapolated ones: 
indeed, on the finite system the gap does not vanish even 
close to the critical locus and consequently the generation 
of defects is smaller than the expected one in the thermo- 
dynamical limit. At small v the extrapolated data are in 
perfect agreement with the scaling prediction Q which 
is represented by the full line. The saturation at large 
v of the defects production, independently of the quench 
protocol (a and uj values), is due to the fact that for very 
fast quenches the only relevant parameters are the initial 
and final amplitudes g. If the initial amplitude is very 
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FIG. 2: Density of defects n versus the quench parameter v for a critical quench. The amplitude changes from g — 1 to g — 0. 
Empty symbols correspond to different system sizes (L = 96 to L — 256 from bottom to top). The extrapolated data (filled 



circles) show, in the adiabatic limit d< 1, a perfect agreement with the scaling prediction: n ~ « 1//4 for to 
for to — 2, a = 2; n ~ v 1 ' 6 for uj = 2, a - 
evaluated on a system with L — 1024 spins. 
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3 (straight lines). The dashed lines give the sudden-quench value n sq ~ 0.179 



high in modulus, the correlation length is very small (of 
order one) and there is no big differences for different 
values of uj. One expects the same defect production as 
in the case of a sudden quench to the critical point of a 
completely disordered initial state (which is given by the 
horizontal dashed line). A similar behavior was reported 
in [13] where for fast enough inhomogeneous quenches 
one recovers the homogeneous defect production (corre- 
sponding here to the homogeneous sudden quench satu- 
ration at large v) while for slow enough inhomogeneous 
quenches (here small v) the defect production is signifi- 
cantly lowered. 

Conclusion We have obtained the scaling properties 
of the density of defects generated during the loading or 
unloading of a power-law trap. The analysis is based on 
the identification of a KZ-time scale which depends on 
the universal properties of the critical point as well as on 
the exponents characterizing the temporal ramp and the 
spatial trap. We found that the optimal nonlinear pas- 
sage through the quantum critical point, specified by an 
optimal time ramping exponent a opt , is strongly affected 
by the presence of the trapping potential. The analyses 
made on the quantum Ising chain revealed quite strong 
finite size corrections, as seen on figure|2j to the expected 
scaling prediction Q for the density of defects. A rele- 
vant extension of this work is the study of the influence of 
a finite temperature on the scaling properties [15] . This 
is currently under investigation. 
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